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We analyze the ghost issue in the recently proposed models of non-linear massive gravity in the 
ADM formalism. We show that, in the entire two-parameter family of actions, the Hamiltonian 
constraint is maintained at the complete non-linear level and we argue for the existence of a non- 
trivial secondary constraint. This implies the absence of the pathological Boulware-Deser ghost to 
all orders. To our knowledge, this is the first demonstration of the existence of a consistent theory 
of massive gravity at the complete non-linear level, in four dimensions. 



(N 

o 

< 



> 

m 
m 

o 



X 



Introduction and Summary: The search for a con- 
sistent theory of massive gravity is motivated by both 
theoretical and observational considerations. Since the 
construction of a linear theory of massive gravity by 
Fierz and Pauli in 1939 [HQ, a proof of the existence of 
a consistent non-linear generalization has remained elu- 
sive, making it a theoretically intriguing problem. On 
the observational side, the recent discovery of dark en- 
ergy and the associated cosmological constant problem 
has prompted investigations into long distance modifica- 
tions of general relativity. An obvious such modification 
is massive gravity. 

Theories of massive gravity generically suffer from a 
ghost instability. The origin of this problem can be un- 
derstood as follows. In general relativity the four con- 
straint equations of the theory along with four general 
coordinate transformations remove four of the six prop- 
agating modes of the metric, where a propagating mode 
refers to a pair of conjugate variables. The total number 
of propagating modes is thereby reduced to the physi- 
cal two modes of the massless graviton. In contrast, in 
massive gravity the four constraint equations generically 
remove the four non-propagating components of the met- 
ric while the general covariance is broken. Thus the the- 
ory will contain six propagating modes of which only five 
correspond to the physical polarizations of the massive 
graviton. The remaining mode is a ghost. 

The question then is whether it is possible to construct 
a theory of massive gravity in which one of the constraint 
equations and an associated secondary constraint elimi- 
nate the propagating ghost mode instead. The linear 
Fierz-Pauli theory succeeds in eliminating the ghost in 
this way. But Boulware and Deser 0] showed that the 
ghost generically reappears at the non-linear level. More 
recently progress was made in Q by observing that the 
ghost is related to the longitudinal mode of the Gold- 
stone bosons associated with the broken general covari- 
ance. This greatly simplifies the analysis of the ghost 
problem in the so-called decoupling limit which isolates 
non-linear effects in the ghost sector. Based on this ap- 
proach a procedure was outlined in 0, HJ to avoid the 
ghost order-by-order by tuning the coefficients in an ex- 
pansion of the mass term in powers of the metric per- 



turbation and of the Goldstone mode. In 2010 dc Rham 
and Gabadadze Q successfully obtained such an expan- 
sion which is ghost free in the decoupling limit. Later in 
0], these perturbative actions were resummed into fully 
non-linear actions resulting in a two-parameter family of 
theories. This was the first successful construction of po- 
tentially ghost free non-linear actions of massive gravity. 
Also in |7| , one of these resummed actions was analyzed 
in the ADM formalism Q and it was argued to be ghost 
free to fourth order in metric perturbations around flat 
space. In Q it is claimed that the ghost still appears at 
the fourth order. (For a review of recent developments 
in massive gravity see fioj].) The present work addresses 
the ghost issue at the non-perturbative level. 

The systematics and generality of these potentially 
ghost free actions are studied in [11|. In particular, they 
are presented as a two-parameter generalization of a min- 
imal extension of the Fierz-Pauli theory. In this work we 
show that the entire two-parameter family of actions is 
ghost free at the full non-linear level. Starting with the 
minimal theory in the ADM formalism, we show that 
the lapse TV is indeed a Lagrange multiplier leading to a 
Hamiltonian constraint on the propagating modes. We 
also show that the same analysis extends to the full two- 
parameter generalization of the minimal theory. We then 
argue that this Hamiltonian constraint gives rise to a sec- 
ondary constraint. These are enough to eliminate a sin- 
gle propagating mode, ensuring that the theory contains 
only five propagating degrees of freedom appropriate for 
the spin-2 massive graviton. Thus the Boulware-Deser 
ghost is eliminated. 

Non-Linear Massive Gravity: In the Fierz-Pauli 
theory, linearized general relativity in fiat space is ex- 
tended by the addition of a mass term for the metric 
fluctuations h^„ = g^ u - rj^ u , 



(1) 



To construct non-linear generalizations of the Fierz-Pauli 
mass term, an additional non-dynamical metric is in- 
variably required. In the recently developed potentially 
ghost free theories, the basic building block is a matrix 
of the form \/g~ l f @, EH , where the square root of the 
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matrix is defined such that \J g~ l f \J g^ 1 f = g^ x f\v In 
particular, the minimal extension of the Fierz-Pauli ac- 
tion with zero cosmological constant is given by fill ]. 



S 



R-2m 2 (trVfl _1 /-3) 



(2) 



Our ghost analysis is based on this action for the case 
of a flat fp V so that in the physical gauge / M „ = ?/ M „. 
Recent studies of massive gravity have primarily focused 
on this case [1, 0] ■ 

However, with this choice of / M „, the minimal action 
does not have a Vainshtein mechanism 12| and thus ex- 
hibits the vDVZ discontinuity [lH, [H| , as shown in [lj| . 
To be compatible with observations one must consider 
theories with additional higher order interactions that 
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could induce a Vainshtein mechanism (see, e.g. 
We will show that our analysis naturally extends to the 
entire family of such actions without any modification. 

The most general non-linear massive gravity theo- 
ries that are potentially ghost-free are given by a two- 
parameter family of actions. Defining a matrix IK so that 
\J g~ l f = 1 + IK, these actions can be written as 0, 



S = Mi 



d A x, 



R+2m 



2 E< 

71=2 



(3) 



with oil = 1 and where the are defined in (|6]) below. 

For our purposes it will be easier to work with an equiv- 
alent formulation of ([3]) 11| . 



S = Ml / d 4 x 



i? + 2m 2 ^/3„ e , 1 (Vg rT 7) 



, (4) 



where the fi n are given in terms of the a n of ([3]) as, 

Po = 6 — 4a 3 + a 4 , f}\ = — 3 + 3a 3 — a 4 , 

fa = 1 - 2a 3 + a 4 , f3 3 = a 3 - a 4 ■ (5) 

The efc(X) are elementary symmetric polynomials of the 
eigenvalues of X. For a generic 4x4 matrix they are 
given by, 



e 2 (x) = i(w - n 

7 3i 



eo(X) = l, ex(X) = [) 
e 3 (X) = ±([X] 3 ~3[X][X 2 ]+2[X 3 ]) 
e 4 (X) = Jj([X] 4 - 6[X] 2 [X 2 ] + 3[X 2 ] 2 + 8[X][X 3 ] - 6[X 4 ]) , 
e fe (X) = for k > 4, (6) 

where the square brackets denote the trace. The action 
(UJ contains terms that are at most third order in \J g~ x f 
rather than fourth order as in ([3]). When a 3 = a 4 = one 
obtains the resummed theory for which the ghost analysis 
was performed in [tJ to fourth order. When a 3 = a 4 = 1 
one obtains the minimal action ([2]). After treating the 
minimal action, we will extend the ghost analysis to the 
most general case (01, for arbitrary (3 n . 



The Hamiltonian constraint: Let us recapitulate 
the counting of degrees of freedom in standard masslcss 
general relativity. In the ADM formulation Q, the ten 
components of the metric are parametrized as 



N= (-g°T 1/2 , 



N, 



goi 



Jij — gij ■ 



(7) 



The jij describe six potentially propagating modes. The 
action written in terms of canonical variables is linear in 
the non-propagating modes N and iVj (collectively, N^). 
Thus the equations of motion are constraints on the 
7y and their conjugate momenta n 11 . Along with the 
general coordinate transformations they eliminate four 
out of six propagating modes, a propagating mode re- 
ferring to a pair of conjugate variables. The are de- 
termined by the remaining equations, thus leaving two 
propagating modes corresponding to a spin-2 graviton. 

In a generic non-linear extension of massive gravity, the 
mass term depends non-linearly (but still algebraically) 
on the N^. The corresponding equations of motion deter- 
mine these non-dynamical variables in terms of jij and 
7r y , keeping all six of the propagating modes undeter- 
mined. Five propagating modes of jij correspond to the 
massive graviton, while the sixth one is a ghost, called 
the Boulware-Deser mode A ghost free theory of 
massive gravity must maintain a single constraint on 7^ 
and 7r lJ along with an associated secondary constraint to 
eliminate this ghost-like sixth mode. Below we show that 
this is indeed the case for the non-linear massive gravity 
actions described above. 

Let us first consider the minimal massive gravity action 
((2|). In the ADM parameterization the Lagrangian C is 
given by, 

i: lJ d tllJ +NB +N l R l -2m 2 ^N (tr - 3) , (8) 

where (with N' 1 = j ij Nj), 

r -1 V* — f 1 Nl5 *i \ ^ 

{9 m v N 2 \-N l (iVV - ATW')% J ' { ) 



Here and in what follows we use y^y to denote yj det 7^ . 

The action © is highly nonlinear in and thus it 
might appear that there are no constraint equations for 
the propagating degrees of freedom. However, if the four 
Nfj, equations only depend on three combinations of N 
and Ni, the fourth equation can be used to determine 
the sixth mode of 7^ in terms of remaining modes. 

To show that is the case, we start by assuming that 
three such combinations n l exist. Then, after writing N l 
in terms of n l , the massive gravity actions should satisfy 
the following two criteria, 

1 . The action is linear in N so that the N equation of 
motion becomes a constraint on the other fields. 

2. The equations of motion for the n l arc independent 
of N and hence are algebraically solvable for the n l . 
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Thus the A equation becomes a constraint on the 7^ and 
Tt lJ . Along with a secondary constraint, this removes the 
ghost. A itself is non-dynamical and is expected to be 
determined by the remaining equations, as in GR We 
will see that when criterion 1 is satisfied, 2 will follow 
automatically. 

Criterion 1 means that the change of variables must be 
linear in A, hence we consider, 

N z = {5) + ND i j )n j . (10) 

The matrix D = D* is determined by requiring that the 
mass term is linear in N. Indeed that will be the case if 
the square-root matrix has the form, 



The symmetry property of D along with expression f| 16[) 
can now be used to show that the n l equations of mo- 
tion are independent of A as demanded by criterion 2. 

Indeed, using 8 tr y/M = | tr(y/M SM) to differentiate 
the trace term, one gets, after some manipulations, the 
n k equation of motion, 



Ri 



i/l — n r 8 rs n s 



= 0. 



The expression in the square brackets is the Jacobian of 
the transformation ([T0[) and is non-zero. Hence the n % 
equations are, 



Nyfg-h) = A + AB , (11) 
where matrices A and IB are independent of A. Then, 

- B 2 . (12) 



(y/l -n r 8 rs n s ) Ri + 2m 2 ^n l 8 u = . 



(19) 



.9 S = 



1 A 2 1 



iV 



These can be readily solved to determine n % in terms of 
jij and the conjugate momenta n 13 , 



N 2 



-RjS ji [4m 4 det 7 + R k 5 kl Ri] 



-1/2 



(20) 



On the other hand, to write g~ x r\ in terms of the new 
variables, let's assemble the n % into a column vector n, 
with transpose n T , and write 77 = diag{— 1, 1}, where, 



This solution implies that y/l — n T In is real. 
The A equation of motion is, 



I = <% , I 1 = 5 ij , whereas 1 = 5). (13) R° + R^)^ - 2m 2 ^7 \y/l- n r S rs n s D 



Then, writing (J9j> in terms of the variables (|10[) and iden- 
tifying the resulting expression with p^|) . we read off, 



1 



1 



y/l - n T In \ ~ n 



-nn T I 



^(7-! - Dnn T D T )I J ' 



(14) 



(15) 



( y/l - n T I n) D = , (j- 1 - Dnn T D T )l . (16) 



and 



This last equation can be easily solved for £)*•. However, 
for the arguments that follow we only need the equality 
(fT6"| and not the explicit solution. Note that the trans- 
formation (|10[) contains no time derivatives and can be 
shown to be invertible. 

A crucial property of D is that is symmetric, 



= 0. 
(21) 

Using the n l solution, this clearly becomes a constraint 
on the 12 components of 7^ and it 13 . Note that in the 
limit that m 2 -> 0, (fl9|) and ([21]) reduce to the four 
constraints of general relativity. 

The general action: We now extend the analysis of 
the previous section to the full two-parameter general- 
ization of the minimal theory. First consider the next 
higher term in \/ g^ 1 ^ in the action given by, 

l r, 



e 2 { 



W 9 1 V) = o ( tr Vg 1 V) 2 -trg 1 T] 



(22) 



To express this in terms of the variables defined in the 
previous section, note that the matrix A has the property 
tr(A fc ) = (tr A) fc . The potential ([22} then gives, 



1 , , 

Ne 2 = - [2(tr Atrl 



tr AB)+A^ ((tr B) 2 - trB 2 )] 



-l\T 



DI~ L = pi" 1 ) 



(17) 



This can be seen from ([T5[) combined with the relation 
(y/MI)!- 1 = \- l (y/lM) which holds for any matrix M ij . 

Now in terms of the new variables n z , the action is 
linear in A, meeting the first criterion, 

C =^d tll] + NR° + RM + ND l An j 



(23) 

This is linear in the lapse N and thus also satisfies our 
first criterion. Varying with respect to n k gives, 



- (rrot 



(Ne 2 ) 



l SuD r 



l 6 lm D™) 



(24) 



2m 2 ^/7 y/l — n T Ir 



-Atr(v / 7" 1 I - Dnn T D T l) - 3A 



(18) 



It is straightforward to show that the next term in the 
potential, A e^{yj g~ 1 vi) 1 is also linear in the lapse A, 
and, through a more involved analysis, determine the 
corresponding contribution to the n k equation of motion. 
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Combining these results with those from the previous 
sections, the complete equations of motion for n l are, 



V v 1 — n r o rs n s 
+/3 3 (Vl - n r <5 rs n s ) n% fc 



D k D m 

i rn 



0. 



(25) 



As in the previous sections, these equations are inde- 
pendent of N and can be used to eliminate n\ The N 
equation is then the constraint on 7^ and ir lJ , 



+ 2m 2 ^[/3 + ft tr B + \fa {(tr B) 2 
+ i/3 3 {(tr B) 3 -3trBtr B 2 + 2 tr B 3 } 



(26) 



tr 



= 0. 



The secondary constraint: We now argue that the 
Hamiltonian constraint gives rise to a secondary con- 
straint (for a proof, see [la ], completed while this work 
was in review). This implies that the 12 dimensional 
phase space of the dynamical variables 7.^ and n 11 has 
only 10 degrees of freedom, corresponding to the five po- 
larizations of the massive graviton. 

We have shown that, upon integrating out the shift 
N*, the Lagrangian (|18p remains linear in the lapse N, 



non-trivial secondary constraint at the non-linear level 
as well. Moreover, as can be seen from the Fierz-Pauli 
structure, enforcing |C, 2) , H } w will result in an equa- 
tion for N, rather than a tertiary constraint. Thus no 
further degrees of freedom are removed in this way. 

Discussion: This work demonstrates the existence of 
non-linear theories of massive gravity that do not suffer 
from the Boulware-Deser ghost instability. Note that, in 
order not to violate the constraints found above, the cou- 
pling of the metric to matter must also be linear in the 
lapse and shift functions. The minimal coupling of Gen- 
eral Relativity automatically satisfies this requirement 
and hence will not change the arguments presented here. 

It should be emphasized that while it is common to 
discuss the ghost in terms of Stiickelberg fields 0-0] i the 
Boulware-Deser instability Q is, strictly speaking, due to 
the loss of the Hamiltonian constraint. We have shown 
that the massive actions ([J) precisely avoid this problem. 
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C = tt« - ^0 (7« , ^ 3 ) + N C ( 7l , , n ij ) • (27) 

A secondary constraint is obtained by demanding that 
the primary constraint C is independent of time on 
the constraint surface. In the Hamiltonian formula- 
tion this condition is given in terms of the Poisson 
bracket, {C, H} « 0, where H = J d 3 x(H - NC). If 
{C(x), C(y)} w 0, then this condition is independent of 
N and thus becomes a constraint on 7^ and tt 1,3 , 



C m = {C, H o }^0. 



(28) 



where now Hq = J d 3 xT-Lo. 

By construction, the Lagrangian (|27p reproduces the 
Fierz-Pauli Lagrangian at lowest order in the fields, 

H^U^ + O^y), C~C^ + 0(7 2 ,7r 2 ). (29) 



Hence one can compute, 



c. 



:C£f + 0(7V 2 ), 



(30) 



where is neither identically zero nor equal to 

C FP . Now, in the Fierz-Pauli case, we know that 
{C FP (x), C FP (y)} w 0. Thus at lowest order in the 
fields there exists a non-trivial secondary constraint. As 
long as {C(x), C(y)} w continues to hold at the non- 
linear level (for a proof, see [IH): then C (2) remains a 
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